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Abstract. We present iV-body simulations of Low Sur- 
face Brightness (LSB) galaxies and their Interstellar Medium 
to investigate the cause for their low star formation rates 
(SFR). 

Due to their massive halos, stellar disks of LSB galaxies 
are very stable and thin. Lack of dust makes the projected 
edge-on surface brightness of LSB galaxies comparable to 
the projected edge-on surface brightness of dust-rich High 
Surface Brightness (HSB) galaxies of similar size. 

We show that the low surface densities found in LSB 
galaxies are by themselves not enough to explain the slow 
evolution of LSB galaxies. A low metal content of the gas 
is essential. As a consequence the gas cools inefficiently, 
resulting in an almost negligible cold gas fraction. We show 
that LSB galaxies must have molecular gas fractions of less 
than 5 percent. 

Our best model has a SFR which is on average low but 
fluctuates strongly. This causes the large spread in colors 
of LSB galaxies. From a distribution of birthrate param- 
eters we conclude that the presently-known and modeled 
gas-rich blue LSB galaxies constitute the majority of the 
total population of gas-rich LSB disk galaxies. We deduce 
the existence of an additional red, quiescent and gas-rich 
population which constitutes less than 20 percent of the 
total population. This does not rule out the existence of a 
large number of gas-poor LSB galaxies. These must how- 
ever have had an evolutionary history dramatically differ- 
ent from that of the gas-rich galaxies. 
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drodynamics — methods: numerical 



1. Introduction 

Deep surveys of the night-sky have uncovered a large pop- 
ulation of disk galaxies with properties quite different from 

* Present address: School of Physics, University of Mel- 
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those of the extensively studied "normal" high surface 
brightness (HSB) galaxies. These so-called Low Surface 
Brightness (LSB) galaxies have, as their name already 
implies, surface brightnesses much lower than what was 
until quite recently assumed to be representative for disk 
galaxies (Freeman 1970). 

The LSB galaxies we will be discussing are generally 
dominated by an exponential disk, with scale lengths of 
a few kpc. Morphologically they form an extension of the 
Hubble sequence towards very late- type galaxies. Observa- 
tions suggest that LSB galaxies are unevolved galaxies, as 
shown by e.g. the sub-solar metallicities (McGaugh 1994) 
and the colors (McGaugh & Bothun 1994; de Blok et al. 
1995; van den Hoek et al. 1997). 

The evolutionary rate of a galaxy may in fact be re- 
flected in its surface brightness. For example, the gas frac- 
tion (Mgas/Mgas+stars) increases systematically with sur- 
face brightness, from a few percent for early type spirals 
to much higher values approaching unity for late type LSB 
galaxies ( |de Blok fc McGaugh 199^ ). In many LSB galax- 
ies the gas mass exceeds the stellar mass (even for such 
extreme assumptions as maximum disk). 

It is still unclear what the physical driver is for the 
difference between HSB and LSB galaxies. Investigations 
of their dynamics, using Hi observations (de Blok et al. 
1996), suggest that LSB galaxies are low-density galaxies 
(de Blok & McGaugh 1996). This is one of the favored 
explanations for the low evolution rate of LSB galaxies 
(see e.g. van der Hulst et al. 1987), as this implies a large 
dynamical time-scale. Environment may also play a role. 
Tidal interactions with other galaxies increase star forma- 
tion rates. This may not happen in LSB galaxies which 
are found to be isolated (Mo et al. 1994). In LSB galaxies 
a few star forming regions are usually found. These are 
distributed randomly over the galaxies and do not trace 
the spiral arms. 

The main purpose of this paper is to investigate whether 
the low density alone can explain the properties of LSB 
galaxies. We use 3D numerical simulations to address this 
problem. The dynamical simulations include both stars 
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and gas and incorporate a parameterized description of 
star formation, including feedback on the gas. 

The prescription of star formation in numerical simu- 
lations is not straightforward, given the limited resolution 
of the models and the fragmentary knowledge about the 
physics governing star formation on kpc scales. In the lit- 
erature various different star formation algorithms can be 
found. For instance, Friedli & Benz (1995) use a crite- 
rion based on the Toomre stability parameter Q (Toomre 
1964) to study star formation in barred systems; Mihos 
& Hernquist (1994a, b, 1996) adopt a Schmidt law based 
on the gas density [star formation rate (SFR) oc p 1 - 5 } to 
study mergers of galaxies; others require gas to be in a con- 
vergent and Jeans unstable flow to form stars flKatz 1992 ; 
Katz et al. 1996[ [Navarro fc White 1993|; jBtcinmctz 1996). 



We employ the method of Gerritsen & Icke (1997, 
1998), which uses a Jeans criterion to define star forming 
regions, coupled with an estimate of the cloud collapse 
time. In these simulations gas is treated fully radiative 
with allowed temperatures between 10 K and 10 7 K; cool- 
ing is described by standard cooling functions, heating is 
assumed to be provided by far-ultraviolet (FUV) radiation 
and mechanical heating from stars. In this way the sim- 
ulated interstellar medium (ISM) mimics a multi-phase 
ISM ( |Ficld ct al. 1969fc |McKee fc Ostriker 1977| ). This is 
an improvement over other simulations which try to cre- 
ate a multi-phase ISM but do not allow radiative cooling 
below 10 4 K (e.g. [Hernquist 1989| ; |Katz et al. 1996; ). The 
multi-phase ISM allows us to restrict ourselves to consid- 
ering only relatively cold (T < 10 3 K) regions as the sites 
for star formation. 

Here we apply this method to study the ISM and star 
forming properties of LSB galaxies, and our results will 
be valid under the assumption that the physical processes 
governing star formation are the same for HSB and LSB 
galaxies. We test whether a low mass-density is sufficient 
to explain the properties of LSB galaxies by discussing two 
different models of a specific LSB galaxy: the first model 
has the cooling properties of a solar abundance ISM, while 
in the second model the cooling efficiency is lowered, thus 
mimicking a metal-poor ISM. 

The structure of the paper is as follows. In Sect. || we 
describe the numerical techniques. The construction of a 
model LSB galaxy is presented in Sect. |[ This galaxy 
model is evolved for a few Gyr, using two different pre- 
scriptions for the cooling properties of the gas (Sect. |J). 
The implications of the simulations are discussed in Sect. 
|5|. We conclude this paper with a summary in Sect. §. 



2. Numerical technique 

We model the evolution of galaxies using a hybrid N- 
body/hydrodynamics code (TREESPH; Hernquist & Katz 
1989). An extensive description is given in Gerritsen & 
Icke (1997, 1998). A tree algorithm (|Barncs fc Hut 1986; 



Hernquist 1987) determines the gravitational forces on the 
collision-less and gaseous components of the galaxies. The 
hydrodynamic properties of the gas are modeled using 
smoothed particle hydrodyn amics (SPH) (see Lucy 1977 ; 
|Gingold fc Monaghan 1977 ). The gas evolves according to 
hydrodynamic conservation laws, including an artificial 
viscosity for an accurate treatment of shocks. Each parti- 
cle is assigned an individual smoothing length, h, which 
determines the local resolution and an individual time 
step. Estimates of the gas properties are found by smooth- 
ing over 32 neighbors within 2h. We adopt an equation of 
state for an ideal gas. 

We allow radiative cooling of the gas according to the 
cooling function for a standard hydrogen gas mix with a 
helium mass fraction of 0.25 (Dalgarno & McCray 1972). 
Radiative heating is modeled as photo-electric heating of 
small grains and PAHs by the FUV field ( |Wolfire et al. 1995| ) 
produced by the stellar distribution. 

All our simulations are advanced in time steps of 3x 10 6 
yr for star particles, while the time steps for SPH particles 
can be 8 times shorter (3.8 x 10 5 yr). The gravitational 
softening length is 150 pc. A tolerance parameter 9 = 0.6 
is used for the force calculation, which includes quadrupole 
moments. 



2.1. Star formation and feedback 

Gerritsen & Icke (1997, 1998) extensively describe the 
recipe for transforming gas into stars and the method for 
supplying feedback onto the gas. The recipe works well for 
normal HSB galaxies, with the energy budget of the ISM 
as prime driver for the star formation. The simulations 
allow for a multi-phase ISM with temperature between 
10 < T < 10 7 K. This allows us to consider cold T < 10 3 
K regions as places for star formation (Giant Molecular 
Clouds in real life). Below we summarize the important 
ingredients of the method. 

From our SPH particle distribution we select conglom- 
erates where the Jeans mass is below the mass of a typical 
Giant Molecular Cloud. In the simulations performed here 
we use M c = 10 6 Mq for this aggregate mass (the method 
is not very sensitive to the exact value for this mass). The 
Jeans mass depends on the local gas properties and is cal- 
culated from the SPH estimates of the density p and sound 
speed s, 

with G the constant of gravity. During the simulations the 
maximum number density achieved is of the order 1 cm~ 3 . 
It follows that only regions below 10 3 K are unstable and 
may form stars. 

We follow unstable regions during their dynamical and 
thermal evolution and if an SPH particle resides in such a 
region longer than the collapse time, 

t u >t c = - JU=, (2) 
^AtyGp 
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half of its mass is converted into a star particle. Exper- 
iments with a different number of particles and different 
star formatin efficiencies show no dependence on these pa- 
rameters, as already shown in Gerritse & fcke (1997). 

Important for our calculations is that we consider star 
particles as stellar clusters with an age. Thus for each 
individual star particle we can attribute quantities like 
the SN-rate, the mass loss, and the FUV-ffux, according 
to its age. We use the spectral synthesis models of Bruzual 
& Chariot (1993) to determine these quantities, where we 
adopt a Salpeter Initial Mass Function (IMF) with slope 
1.35 and with lower and upper mass limits of 0.1 Mq and 
125 Mq respectively. 

The radiative heating for a gas particle is calculated by 
adding the FUV-flux contributions from all stars, which is 
done together with the force calculation. The mechanical 
luminosity from a star particle is determined by both the 
SN-rate and the mass loss rate. We assume that each SN 
injects 10 51 ergs of energy and that the energy injected by 
stellar winds is E w — ^Am^v^, with the wind termi- 
nal velocity and Am* the stellar mass loss. For massive 
stars Voo depends critically on the stellar mass, luminosity, 
effecti ve temperature, and mcta llicity (e.g. Leitherer et al. 
1992; Lamers fc Leitherer 1993 ). For simplicity we adopt 
Voo = 2000 km/s. After 3.3 x 10 7 yr the last 8 Mq stars 
explode and no more mechanical energy is supplied to the 
gas. Thus we return mechanical energy from massive stars 
into the ISM, and ignore the mechanical energy from low 
mass stars. 

In the simulations the parent SPH particle is the car- 
rier of the mechanical energy from the new star particle. 
This SPH particle ( "SN particle" ) mimics a hot bubble in- 
terior. Radiative cooling is temporarily switched off (the 
resolution does not allow the creation of a low-density, hot 
bubble), the temperature of the SN particle is set to the 
mechanical energy density (of a few 10 6 K), and the par- 
ticle evolves adiabatically. The first 10 7 yr, the position 
and velocity of the SN particle are equal to the position 
and velocity of the associated star particle. Afterwards, 
the SN particle evolves freely. After 3 x 10 7 yr radiative 
cooling is switched on again, and the particle behaves like 
an ordinary SPH particle. 



3. Modeling an LSB galaxy 

There are many ways to construct model galaxies. For our 
purpose we prefer to build a galaxy model after an existing 
galaxy with well-determined properties. Here we choose 
LSB galaxy F563-1; this galaxy is a late- type LSB galaxy, 
representative of the field LSB galaxies found in the sur- 
vey by Schombert et al. (1992). The optical properties of 
this galaxy are described in de Blok et al. (1995); mea- 
surements of metallicities in H n regions are described in 
dc Blok & van der Hulst (1998a); a neutral hydrogen map 
and rotation curve are given in de Blok et al. (1996). Pa- 
rameters such as stellar velocity dispersion, which cannot 



Table 1. Parameters for F563-1 (H = 75 km s -1 Mpc" 1 ). 
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be measured directly, are set in comparison with values 
measured locally in the Galaxy. The current star forma- 
tion rate as deduced from Ha imaging is taken from van 
den Hoek et al. (1997). For convenience these data are 
summarized in Table ^. 

The most difficult problem we face in constructing a 
model is converting the measured luminosity to a stel- 
lar disk mass. This is one of the most persistent prob- 
lems in analyzing the dynamics of galaxies, and, unfortu- 
nately, the present observations do not provide a unique 
answer for this stellar disk mass-to- light ratio [M/Lb]*- 
Rather than using the so-called "maximum disk" value 
(M/Lb)* — 9, which is an upper limit to the possible 
values of (M/Lb)*, we adopt a value based on colors and 
velocity dispersions of (M/ Lb)* = 1.75. An extensive mo- 
tivation for this choice is given by de Blok & McGaugh 
(1997). The implications of choosing a different value of 
(M /Lb)* for the evolution of the stellar disk will be dis- 



cussed in Sect. 3.2 



3.1. The model 

The stellar disk particles are distributed radially accord- 
ing to the (measured) surface density profile. We adopt a 
vertical distribution of the form sech 2 (z/z») (e.g. van der 
Kruit & Searle 1982), with a constant vertical scale height 
z* = 0.5 kpc The disk is truncated at 25 kpc. The lumi- 
nosity yields a total stellar mass of M* — 2.36 x 10 9 Mq. 
The velocity dispersion of the stars is fixed via the relation 

a z = V7rGz*XU, (3) 

where £* is the stellar surface density. This implies a sta- 
bility parameter 3 < Q < 3.5 throughout for the disk 
(dotted line in Fig. |l|). 

Particle velocities are assigned according to the (gas) 
rotation curve corrected for asymmetric drift. Dispersions 
in R,z,6 directions are drawn from Gaussian distributions 
with dispersions of an,a z , ag respectively, using the rela- 
tions valid for the solar neighborhood 



0.6ctr, o~q — —o-r, 

2uj 



(4) 



with u> and k the orbital and epicyclic frequencies respec- 
tively. 
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The gas surface density distribution is modeled after 
the Hi distribution. The total gas mass is M g = 3.85 x 
10 9 Mq (this is the total Hi mass multiplied by 1.4 to 
include He); the surface density decreases almost linearly 
with radius out to 33 kpc. The vertical distribution is as- 
sumed to decline exponentially. The scale height z g of the 
gas layer can be calculated using 

a 2 = 2nG(Z g + E* 



2 -2* ~\~ 



(5) 



(Dopita & Ryder 1994) where a g = 3 km/sis the (adopted) 



gas velocity dispersion. The gas particles are assigned ve- 
locities according to the rotation curve, with isotropic dis- 
persion <Jg. 

The halo is included in the calculations as a rigid po- 
tential. This is justified since the galaxy model evolves in 
isolation. We will thus also ignore any contraction of the 
halo under the influence of the disk potential. We do not 
expect this effect to be important anyway as the mass of 
the disk (assuming (M/Lb)* — 1.75) is only 4 per cent 
of the measured halo mass (de Blok & McGaugh 1997). 
As advantages we do not have to make assumptions about 
halo particle orbits, and we do not have to spend time in 
calculating the force of the halo particles on the galaxy. 
We assign an isothermal density distribution to the halo, 
Po 



Ph 1 + (r/R c y ' 
with central volume density po 



(6) 



0.0751Mq/pc 3 and core 
radius R c — 1.776 kpc. These are the values derived from a 
rotation curve decomposition assuming (M/Lb)* — 1.75. 
This will correctly put the maximum rotation velocity at 
113 km/s. 

For the simulations we use 40,000 SPH particles and 
80,000 star particles initially. This corresponds to an SPH 
particle mass of 9.6 x 10 4 Mq and a star particle mass of 
3.0 x 10 4 Mq. 

3.2. Implications of (M/Lb)* for star formation history 

We can show that a maximum disk value for (M/Lb)* 
is not a plausible option for our modeling exercise, and 
that in fact the "most likely" value we use in constructing 
models may still overestimate the true value of (M/Lb)*, 
making the stellar disk a really insignificant component of 
the whole galaxy system. 

The star formation history of galaxies is commonly 
parameterized by an exponential function: 

SFR = -M g = Mg/T* (7) 

where t* is the star formation t ime scale (e.g. Guiderdon i 
& Rocca-Volmerange 1987 and Chariot fc Bruzual 1991 ). 
This ignores any gas locked up in long-lived stars, but this 
effect will not be important in LSB galaxies due to their 
large gas- fractions. A small value for t* means that star 
formation proceeds rapidly, while a large t* indicates that 
star formation proceeds very slowly. 
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Fig. 1. The stability parameter Q versus radius. The dot- 
ted line shows the Q-value upon initialization of the galaxy 
model. The solid line denotes simulation L, the dashed 
line simulation H. Both simulation yield approximately 
the same stability for the stellar disk 



For HSB galaxies the star formation time scale is nor- 
mally a few Gyr (Kennicutt et al. 1994). For F563-1 the 
measured current SFR as derived from Ha imaging is ap- 
proximately 0.05 Mq /yr; this yields a star formation time 
scale of 77 Gyr. Thus star formation in this disk proceeds 
much slower than in HSB galaxies. 

Since we know the stellar and gas mass, and hence the 



total disk mass M t = M* + M g of the galaxy, we can use 
the exponential parameterization and the known SFR to 
compute the "age" td of the disk (i.e. the time elapsed 
since star formation started): 



M* = 



M g dt = M t ( 



1 



-t d /r. 



(8) 



For the adopted stellar mass of F563-1 Eq. || yields an age 
for the galaxy of 37 Gyr, much larger than the age of the 



universe. 



If the exponentially declining SFR describes the star 
formation history adequately and the galaxy is not in a 
state of unusually low star formation activity (note that 
the blue colors for these LSB galaxies suggest a relatively 
high current star formation activity) , then the old implied 
age of the disk could mean that the mass of the disk should 
in reality be even less. 

We show this in Fig. || where we present a few solutions 
for Eqs. ^ & ^ for various values of the SFR. If we adopt 
0.05 M /yr as the true SFR then (M/Lb)* ~ 0.6 if the 
galaxy is 15 Gyr old. Higher values for (M/Lb)* would 
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Fig. 2. Relation between galaxy age, td (vertical axis), and 
stellar mass, M* (horizontal axis), according to Eqs. & 
|8|. The different lines are solutions to these equations for 
various values of the current SFR, where the measured 
value is 0.05 Mq/jv. The top axis shows the (M/Lb)* 
ratio corresponding to the stellar mass. The input model 
assumes (M/Lb)* = 1.75, the maximum disk model has 
(M/Lb)* = 9, which yields a very high age for the galaxy 



make the galaxy older than the universe. A maximum disk 
solution (van Albada & Sancisi 1986) for the stellar mass- 



to-light ratio yields (M/Lb)* = 9, which leads to an age 
of 110 Gyr. A current SFR of order 0.4 Mq /yr is needed 
to reconcile the maximum disk with the exponentially de- 
clining SFR model, which implies that the measurements 
underestimate the true SFR by an order of magnitude. 
Thus seems highly unlikely (it would give LSB galaxies 
SFRs comparable to those of actively star forming late 
type galaxies), and therefore we reject the maximum disk 
as an acceptable solution for the stellar disk mass. 

The adopted (M/Lb)* = 1-75 which we use in the 
galaxy model is only consistent with the age of the uni- 
verse if the true SFR is of order 0.10 Mq /yr, otherwise the 
stellar mass is probably overestimated. However if the stel- 
lar disk is indeed less massive than our galaxy model, then 
any conclusions concerning the disk stability and thickness 
of the disk will be even stronger. 

It might be argued that the exponential star formation 
history is not the right choice for this type of galaxy, and 
that other functional forms of star formation history, in 
combination with the measurements, will yield reasonable 
ages. This is not true, as any other (reasonable) star for- 
mation history will result in even larger ages for the disk. 



Taking for example another extreme functionality - a con- 
stant SFR of 0.05 Mq /yr - yields an age t d = 47 Gyr. The 
only way to get a reasonable age, given the measured lu- 
minosity and assumed (M/Lb)*, is to assume that the 
current SFR in F563-1 (and LSB galaxies in general) has 
been underestimated, which is unlikely. We refer to Mc- 
Gaugh & de Blok (1997) for an extensive discussion on 
the various functionalities of the star formation history. 

We therefore conclude that the disk of F563-1 must 
have a small stellar mass-to-light ratio. 



4. Evolution 

In this section we demonstrate that the low density as 
found in LSB galaxies, by itself is not sufficient to repro- 
duce the low observed SFRs of LSB galaxies. Low metallic- 
ity gas is required to explain the properties of LSB galaxies. 

To show this we construct two model galaxies using the 
structural parameters of F563-1. Model H represents an 
LSB galaxy with a solar metallicity gas. Although we use 
the structural parameters relevant for F563-1, the model 
is in effect a model HSB galaxy, which is "stretched out" 
to give the low (surface) densities found in LSB galaxies. 
This model therefore tests the low-density hypothesis. 

The other model, L, has the same structural parame- 
ters as H, but in addition we lowered the cooling efficiency 
of the gas below 10 4 K by a factor of seven. Cooling below 
10 K is dominated by metals, so lowering the efficiency is 
equivalent to lowering the metallicity by an equal amount. 
Model L thus most closely approximates what is currently 
known observationally about LSB galaxies. In summary 
we adopt the cooling function of a solar abundance gas in 
model H . This cooling function is later changed in model 
L to simulate the effects of low metallicity. 

We assume that the physics regulating the star forma- 
tion is the same for LSB and HSB galaxies, and therefore 
use the star formation recipe described in Section 2 and 
applied to HSB galaxies in Gerritsen & Icke (1997, 1998). 
After initialization we let the two model galaxies evolve 
for 2.3 Gyr. The models rapidly settled into equilibrium, 
any longer time interval would have produced identical 
results. 

In this section we first discuss the evolution of the stel- 
lar disk. The stellar disk evolves rather independently from 
the gas disk and star formation. Hence we effectively ex- 
plore the consequences of the disk/halo decomposition, 
notably on the stability and thickness of the stellar disk. 
Second, we discuss the evolution of the SFR with time for 
both simulations. For a physical interpretation of the dif- 
ference in SFRs we present phase diagrams of the ISM in 
both simulations. 
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4--1- Stellar disk 

Figure [l] shows the stability parameter Q for the stellar 
disk at t = 2 Gyr (that is 2 Gyr after the start of the 
simulation which started at t = 0), where Q is defined as 

K(J R 



3.36GE* 



(9) 



with k the epicycle frequency (Toomre 1964). Q < 1 de- 
notes a (local) instability, while Q > 1 means stability. For 
both simulations Q « 2.5 in the center and rises to Q ~ 5 
in the outer parts of the disk. There is no evolution in the 
stability. Once the system has settled after the start of the 
simulation, the radial behavior of Q as shown is reached. 
The value of course depends on the input parameters, no- 
tably the stellar mass, but it is clear that the stellar disks 
of LSBs are more stable than the stellar disks of HSBs, 



where Q is of order 2 ( Bottema 1993 , van der Hulst et al. 
1993). 

In both simulations the scale height of the stellar disk 
decreases by approximately a factor of 1.3. The initial 
scale height corresponded to the scale height of an Sc 
galaxy (van der Kruit & Searle 1982). The final axial ra- 
tio for the galaxy model is about 15. As particle scattering 
during the simulation tends to increase the thickness of the 
disk, we conclude that the stellar disks of LSB galaxies are 
thinner than those of HSB galaxies. We will return to this 
in more detail in Sect. 5.1. 



4.2. Star formation & ISM 

Figure || shows the SFR versus time. In model H the SFR 
settles after a short adjustment at a value of 0.28 Mq /yr 
and declines slowly with time. The SFR in model L settles 
at 0.08 M /yr. 

For both simulations the SFR varies on time scales of 
a few tens of Myr and the amplitude of these variations 
can exceed 0.1 Mq /yr. The rapid variability in the SFR 
is due to the discrete nature of star formation in our sim- 
ulation. New star particles have a mass of approximately 
5 x 10 4 Mq and thus represent (large) stellar clusters. In 
this way our simulations incorporate the idea that most 
stars in real galaxies form in stellar clusters. For a small 
number of star forming regions (as in LSB galaxies) this 
will in real life also lead to a rapidly fluctuating, but on 
average low SFR. 

In simulation L the strength of the star formation peak 
is large compared to the average SFR. These fluctuations 
in star formation activity will actually dominate the color 
of the galaxy, with large fluctuations giving rise to blue 
colors. The same star formation fluctuations will also be 
present in HSB galaxies, but due to the higher SFR the 
large number of fluctuations will give the impression of a 
high, relatively constant SFR. 

We thus find that it is not so much the absolute value 
of the average SFR which determines the colors, but the 



Lll 

in 



0.40 



0.30 



0.20 



0.10 



0.00 



,A'!j 




0.5 



2.0 



1.0 1.5 
Time [Gyr] 

Fig. 3. The evolution of the SFR versus time. The lower 
line shows the SFR from simulation L, the upper line rep- 
resents simulation H . Data are averaged over 9 Myr 



contrast of any SF fluctuation with respect to the average 
SFR. For the LSB galaxies the large contrast leads to blue 
colors. 

In order to understand the behavior of the star forma- 
tion activity in the simulations we constructed phase dia- 
grams of the ISM. Figure f| shows the temperature versus 
density for all gas particles in both simulations. The top 
panel shows simulation L while the bottom panel shows 
the phase diagram for simulation H . Most particles have 
a temperature of 10 4 K. In simulation H a large frac- 
tion of the particles has a temperature of about 100 K. 
This simulation shows a dominant two-phase structure, 
and resembles the ISM in simulations of HSB galaxies 
( Gerritsen fc Icke 199S ) . Quantitatively, the cold gas frac- 



tion (T < 1000 K) makes up 37% of the total gas mass in 
simulation H and only 4% in simulation L. This reflects 
the cooling properties of the gas: in simulation L seven 
times less heat input is required to keep the gas at 10 4 K 
as in simulation H. 

Consequently simulation L contains virtually only a 
warm, one-phase ISM. The simulated absence of metals 
prevents the ISM from cooling efficiently. As molecular 
gas is only formed in/from the cold component of the 
ISM, the 4% mentioned above is actually an upper limit 
to the amount of molecular gas that could form in such a 
galaxy. The amount of cold molecular gas in sites for star 
formation would thus be negligible. From these models we 
expect the disks of LSB galaxies to contain only negligible 
amounts of molecular gas. This is consistent with obser- 
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Fig. 4. Phase diagrams (temperature versus number den- 
sity) for the simulations. Top panel shows simulation L, 
bottom panel shows simulation H; each dot represents an 
SPH particle. The individual particles at the top of each 
diagram are hot SN particles. Simulation H clearly shows 
a two-phase structure, while in simulation L almost all gas 
is in the warm (T » 10 4 K) phase (85%) 



Fig. 5. Edge-on views of the mass and light distribution 
in LSB galaxy L. The top panel shows a linear grey-scale 
representation of the projected mass surface density. The 
bottom panel shows a grey-scale representation of the pro- 
jected V-band surface brightness. Here, the star particles 
formed during the run-time of the simulation are assigned 
a V^-band flux according to their age following Chariot 
& Bruzual (1993). The contrast has been adjusted to a 
dynamic range of 5 magnitudes. This view resembles the 
so-called "chain galaxies" described by Cowie et al. (1995) 



vations by Schombert et al. (1990) and de Blok & van der 
Hulst (1998b), who find upper limits of less than 10%. 



5. Discussion 

The structural parameters of the modeled disks appear to 
be independent from the star formation activity. In this 
section we will discuss the appearance of the stellar disk 
of LSB galaxies, and the ISM in LSB galaxies. 



5.1. Stellar disk 

LSB galaxies are flatter than normal HSB galaxies (cf 
Fig. 2). From a study of the catalog of edge-on galaxies 
by Karachentzev et al. (1993), Kudrya et al. (1994) con- 
cluded that galaxies become increasingly thinner towards 
later Hubble types, with the very thinnest galaxies hav- 
ing major over minor axis ratios of 20. Our model galaxies 
have axial ratios of 15, that is, almost as flat as the flattest 
galaxies in the Karachentzev et al. catalog. The top panel 
in Fig. ^ shows the edge-on view of mass distribution in 
LSB galaxy L. 

Our model LSB galaxies support claims made by Dal- 
canton & Shectman (1996) that some of the extremely flat 
galaxies ("chain galaxies") detected in medium-redshift 
HST images (Cowie et al. 1995) are simply edge-on LSB 
galaxies. Cowie et al. present HST wide-/ band obser- 
vations of these chain galaxies which they describe as 
"extremely narrow, linear structures ... with superposed 
bright 'knots'." A few of these galaxies were found to lie 
at a redshift z ~ 0.5. As an /-band view of that rcdshift 
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corresponds approximately to a V^-band view locally, we 
have attempted to simulate a U-band observation of model 
L edge-on, by assigning to each stellar particle the /-band 
flux corresponding to its age (Chariot & Bruzual 1993). 
This is shown in the bottom panel of Fig. ||, where the 
dynamic range of the picture has been adjusted to corre- 
spond approximately with the contrast that the observed 
chain galaxies have with respect to the sky-background. 
Comparison with any of the galaxies in Fig. 20 of Cowie 
et al. will show the resemblance. 

How will the edge-on surface brightness distributions 
of LSB galaxies compare with those of HSB galaxies of 
similar scale length? Compared to an equally large HSB 
galaxy, the scale height of an LSB galaxy is smaller. Due 
to the low dust content in LSB galaxies (McGaugh 1994), 
effects of edge-brightening will be stronger in the LSB 
galaxy compared to the dust-rich HSB galaxy. An edge- 
on LSB might have a higher apparent surface brightness 
than an edge-on HSB galaxy of the same size. For exam- 
ple, for HSB galaxies in the Ursa Major cluster, Tully et 
al. (1997) find a differential extinction in the -B-band go- 
ing from face-on to edge-on of 1.8 mag arcsec" 2 . The LSB 
galaxies in their sample are found to be almost transparent 
with a differential extinction of only 0.1 mag arcsec" 2 . The 
almost two magnitudes extra extinction in the HSB can 
thus easily compensate the few magnitudes difference in 
intrinsic surface brightness, making edge-on LSB galaxies 
as bright in surface brightness as edge-on HSB galaxies. 

Hence it is much more difficult to distinguish between 
edge-on LSB and HSB galaxies than between face-on LSB 
and HSB galaxies. A smaller number of star forming re- 
gions, extreme flatness and blue colors will be the only 
easily visible distinguishing characteristics. 

The reason why the stellar disks of LSB galaxies are so 
thin is that they are very stable to local instabilities (Fig. 
[l]). As a consequence, bars and spiral structure, which are 
the most efficient methods for heating a stellar disk (Sell- 
wood & Carlberg 1984), are unlikely to develop sponta- 
neously in LSB galaxies. Thus there is no natural way to 
make the disks thick. 

This explanation is supported by the rarity of barred 
LSB galaxies. In the LSB galaxy catalog by Impey et al. 
(1996), only 4 percent of the galaxies are barred, while the 
frequency of barred galaxies in the RC2 is some 30 percent 
( |Elmcgreen et al. 1990| ). Sellwood & Wilkinson (1993) give 
a fraction of 2/3 for barred HSB galaxies. The stability of 
LSB disks is confirmed in a numerical study of the dy- 
namical stability of these systems by Mihos et al. (1997). 
Whereas we adopted the "most likely" solution for the 
mass-to-light ratio of stellar disk, these authors considered 
the "worst case" scenario of maximum disk. However their 
conclusion is the same as ours: the disks of LSB galaxies 
are extremely stable. 

The lack of truly LSB edge-on galaxies (Schombert et 
al. 1992) could tell us something about the existence of 
very LSB galaxies. As shown above, the presently known 



population of LSB galaxies, when turned edge-on, bright- 
ens to surface brightnesses comparable to those of similar- 
sized HSB galaxies. In principle, the edge-on counterparts 
of the currently known LSB galaxies should thus already 
be in the conventional galaxy catalogs, showing up as 
"streaks on the sky" . As truly LSB edge-on galaxies seem 
to be lacking from the LSB catalogs, this suggests that 
galaxies with face-on surface brightnesses Ho(B) > 25 mag 
arcsec -2 are probably very rare. 



5.2. Blue & red LSB galaxies 

Fluctuations in the SFR as shown in Fig. 3, may very 
well explain why most of the LSB galaxies detected in 
surveys are blue. As McGaugh (1996) argues, the selec- 
tion effects against finding red LSB galaxies on the blue 
sensitive plates on which surveys have been carried out 
are quite severe. It will be very hard to find LSB galaxies 
with colors B — V ~ 1 or redder. So if the absence of red 
LSB galaxies from the current catalogs is caused purely by 
these selection effects, the blue LSB galaxies could simply 
be the bursting tip of a proverbial iceberg. 

In order to explain the colors of the bluest LSB galax- 
ies in the sample of de Blok et al. (1996), van den Hoek 
et al. (1997) had to invoke bursts with a duration of be- 
tween 0.5 and 5 Myr and an amplitude between 1 and 5 
Mq/yt. The fluctuations found in our simulation L have 
an amplitude of 0.15 M©/yr, with a duration of about 20 
Myr. This duration is determined by the lifetime for OB 
stars and the collapse time for molecular clouds, which 
are both of the order of a few 10 7 yr. They are therefore 
slightly milder than the bursts invoked by van den Hoek 
et al. (1997), but keep in mind that the latter bursts were 
used to explain the bluest galaxies. Our models attempt to 
simulate an average LSB galaxy, and it should not come as 
a surprise that the fluctuations we find are milder. Given 
the many assumptions and uncertainties in both modeling 
and observations, it is actually quite encouraging that the 
parameters agree to within an order of magnitude. 

Assuming that blue LSB galaxies are currently under- 
going a period of enhanced star formation, implies that 
there exists a population of red, non-bursting, guiescent 
LSB galaxies. These should then also be metal-poor and 
gas-rich, and share many of the properties of the galaxies 
we have been modeling. 

We can estimate the fraction of red LSB galaxies by 
calculating the distribution of the birthrate parameter b 
for simulation L. The b parameter is the ratio of the 
present SFR over the average past SFR and is a useful 
tool for studying the star formation history of galaxies. 
Birthrate parameters have been determined for a large 
sample of spiral galaxies by Kennicutt et al. (1994). The 
trend is that early type galaxies have small values for b, 
thus most star formation occurred in the past, while late 
type and irregular galaxies have large values for b, often 
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galaxy to be "blue" if its birthrate parameter b exceeds 
the average value of b for a HSB late-type galaxy (see 



1 .5 



Fig. 6. (a) Distribution of the birthrate parameter b, the 
ratio of the current SFR to the average past SFR. Solid 
line shows the b values derived for simulation L, dotted line 
shows the b distribution for a simulation of an Sc galaxy 
(Gerritsen & Icke 1998). On the top are average b values 
for different types of galaxies (Kennicutt et al. 1994). The 
bottom panel shows the cumulative distribution for the 
two simulations. Less than 20% of the b values for the 
LSB galaxy are below b = 0.4, hence we expect at most 
20% of LSB galaxies to be "red" 



exceeding 1, indicating that those galaxy are still actively 
forming stars, and more so than in the past. Here we ap- 
ply this analysis to simulation L in order to estimate the 
fractions of blue and red LSB galaxies. 

In Fig. ^|a we plot the distribution of b values for sim- 
ulation L (solid line), where we have followed the value of 
b over the duration of the simulation in steps of 15 Myr. 
Thus if the SFR peaks in a particular time interval, the 
corresponding value of b will be high. If the SFR is low 
in this time interval b is also low. In total we have 200 6 
values. 

Also shown in Fig. |^a are the b distribution for a simu- 
lation of an HSB Sc galaxy ( Gerritsen fc Icke 1998) , dotted 



line) and the mean values for different galaxy types (from 
Kennicutt et al. 1994). Due to the low average SFR the 



distribution for the LSB simulation is much broader than 
the distribution for the HSB simulation, and the average 
b value is larger. 

Also it is clear that the LSB galaxy has b values larger 
than the average for early type galaxies for most of the 
time. "Classical" LSB galaxies are blue compared to HSB 
galaxies, they thus have an excess of recent star formation 
or equivalently a higher b value. We now define a LSB 



Kennicutt et al. 1994 for relations between birthrate pa- 
rameter and color). Fig. ^a shows that this requires that 
hhSB > (b H SB) ~ 0.4. 

LSB galaxies that do not meet this requirement are 
"red" : non-bursting, but nevertheless still gas-rich. From 
Fig. |^b (which shows the cumulative b distribution) we 
can see that over 80 percent of the fluctuations result in 
blue LSB galaxies. Less than 20 percent of the fluctuations 
therefore results in red LSB galaxies. 

We can estimate colors of the red population using the 
burst models from van den Hoek et al. (1997). For exam- 
ple, a 5 Myr burst with an amplitude of 3 M Q /yr super- 
imposed on a 13 Gyr old population which has undergone 
an exponentially decreasing SFR with a timescale of 4 Gyr 
changes colors and surface brightnesses by A(B — V) ~ 
-0.4, A(R-I) ~ -0.2, AB ~ -1 and Al ~ -0.4 (cf. Ta- 
ble 5 in van den Hoek et al.). Using these values together 
with the measured colors of the bluest LSB galaxies in de 
Blok et al. (1995), yields B - V ~ 1, R - I ~ 0.6 and 
Ho(B) = 24.5 for the red population. 

A recent CCD survey (O'Neill & Bothun 1997) has 
picked up a class of LSB galaxies which have Ho{B) ~ 24 
and B — V ^ 0.8. If some of these galaxies are indeed the 
non-bursting counterparts of the blue LSB galaxies, they 
should be metal-poor and gas-rich, and share many of the 
properties of the galaxies we have been modeling. 

In summary, if the blue colors found in LSB galaxies 
are the result of fluctuations in the star formation rate, 
then this implies that the red gas-rich LSB galaxies con- 
stitute less than 20% of the gas-rich LSB disk galaxies. 
This does not rule out the existence of a population of red, 
gas-poor LSB galaxies. These must however have had an 
evolutionary history quite different from those discussed 
here and possibly have consumed or expelled all their gas 
quite early in their life. 

5.3. ISM 



The SFRs from both simulations are larger than the ob- 
served rate of about 0.05 Mq /yr. This latter value is 
probably correct within a factor of 2. Especially the high 
SFR of simulation H (0.33 Mq /yr) is almost an order 
of magnitude larger than what is needed and it seems 
impossible to reconcile this value with the observations. 
Instead the global SFR is consisten t with the globa l SFR 
for equal- luminosity galaxies (e.g. Kennicutt 1983 ). For 
instance the Sc galaxy NGC 6503 with a maximum ro- 
tation velocity of v = 120 km/s has a total measured 
SFR of 0.4 AWyr, an d a simulated SFR of 0.35 M /yr 
( |Gerritsen fc Icke 1997| , 1998). In general the SFRs of equal- 
luminosity galaxies differ by a factor 10 between LSB 



and HSB galaxy (van den Hoek et al. 1997, and Kenni- 
cutt 1983), despite the copious amounts of gas available 
in the LSB galaxies ( de Blok et al. 199^ ) . 
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SFRl approaches the correct value and the physical 
reason is shown in Fig. ^. The essential information to 
retain from this phase diagram is that we need a differ- 
ent ISM for LSB galaxies, where the bulk of the gas is 
not directly available for star formation, as it is too warm 
(of order 10 4 K). The simulations do not include phase 
transitions from neutral to molecular gas, but as an esti- 
mate for the H2 mass we can consider all star forming gas 
(T ^ 300 K) to be molecular. This gas represents less than 
2% of the total gas mass. In this respect it is interesting 
to note that for a small sample of LSB galaxies CO is not 
detected, yielding an limit of M(H 2 )/M(Hl) < 30 percent 
( |de Blok fc van der Hulst 1998b| , Schombert et al. 1990). 
A few galaxies have upper limits less than 10%. 

As discussed in the previous section the warm ISM 
may be caused by a low metallicity for the gas. Below 10 4 
K, cooling is dominated by heavy elements like C + , Si + , 
Fe + , O. If these elements are rare then it is difficult for 
the gas to radiate its energy away. 

Direct observational support for a low metallicity ISM 
in LSB galaxies comes from oxygen abundance measure- 
ments of Hn regions in LSB galaxies. Those studies yield 
mctallicities of approximately 0.5 times solar metallicity 
( McGaugh 1994 ). Measurements of the oxygen abundance 
in F563-1 (de Blok & van der Hulst 1998a) give an average 
oxygen abundance of 0.15 Zq (compare with the difference 
in metallicity between models H and L). 

A point of concern is the IMF. Our premise is that the 
IMF is universal, however it can differ from our adopted 
Salpeter IMF (see 2.1). Especially a low upper- mass cutoff 
could easily lead to an underestimation of the SFR (e.g. 
an upp er mass cutoff at 30 Mq reduces the measured SFR 
by 2-3; Kcnnicutt 1983j ). However current determinations 
of the IMF in external galaxies point towards a univer- 
sal IMF ( [Kennicutt 198$ |Kcnnicutt ct al. 1994j ). Further- 
more, observations suggest that massive stars are present 
in the Hn regions of LSB galaxies (McGaugh 1994). 



6. Conclusions 

What causes the low evolution rate for LSB galaxies? The 
low density has often been invoked to explain this, since 
the dynamical time scales with l/-»/p. This scenario is 
exactly what is tested in simulation H. The only differ- 
ence with a normal HSB galaxy is the scale length of the 
galaxy, which is for instance three times larger than the 
scale length for the equally bright galaxy NGC 6503. The 
result is striking: adopting "standard" values for the star 
formation process results in a SFR identical to the rates 
of HSB galaxies. 

Thus the low density in itself seems not capable of do- 
ing the job, and we have to rely on a scarcity of heavy 
elements to reproduce a true LSB galaxy. This fits in logi- 
cally with the notion that stars are the producers of these 
elements; the low star formation activity prevents metal 



enrichment of the ISM. It implies that the SFR has been 
low throughout the evolution of LSB galaxies, and that 
these galaxies are "trapped" in their current evolutionary 
state: low density prevents rapid star formation, which 
prevents enrichment of the ISM, which prevents cooling, 
resulting in a warm one-phase ISM. So although the lack 
of metals is directly responsible for the low SFR, the low 
density may ultimately determine the fate of LSB galaxies. 

In summary we reach the following conclusions on the 
physical properties of LSB galaxies. 

1. Dominant halo: The dominance of the dark matter 
halos in LSB galaxies (de Blok & McGaugh 1997) re- 
sults in very thin and stable disks. The major over 
minor axis ratios of the stellar disks are larger than 
15. Lack of dust in LSB galaxies will result in edge- 
on LSB galaxies having apparent surface brightnesses 
equal to those of edge-on HSB galaxies of similar size. 
The lack of truly LSB edge-on galaxies may indicate 
that disk galaxies with ^o(B) > 25 mag arcsec -2 are 
rare. 

2. Low metallicity: Modeling the low SFR found in 
LSB galaxies, demands that the metallicity in the ISM 
must be approximately 0.2 solar, which is consistent 
with observation by McGaugh (1994) and de Blok & 
van der Hulst (1998a). These low mctallicities prevent 
effective cooling so that the ISM in LSB galaxies pri- 
marily consists of warm (~ 10 4 K) gas. The amount 
of cold (molecular) gas is probably less than 5% of the 
Hi mass, supporting claims by Schombert et al. (1990) 
and de Blok & van der Hulst (1998b) who derive small 
molecular gas fractions. 

3. Low average SFR: Due to the fluctuations in the 
SFR in LSB galaxies and their large contrast with the 
average SFR, the spread in colors among LSB galax- 
ies will be larger than among HSB galaxies. From the 
distribution of birthrate parameters in our simulations 
we deduce that, if the currently known blue gas-rich 
LSB galaxies are the most actively star forming LSB 
galaxies, they constitute over 80 percent of the total 
population of gas-rich field LSB disk galaxies. This im- 
plies that there is at most an additional 20 percent of 
quiescent, gas-rich LSB disk galaxies. This does not 
preclude the existence of an additional red, gas-poor 
population. However this population must have an evo- 
lutionary history quite different from that described in 
this work. 
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